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INTRODUCTION 


This paper reviews our understanding of the influence of delamination on the 
fracture sequence and final fracture of graphite-epoxy laminates. The emphasis here 
is on tensile fracture. 

The paper is organized as follows: 

(1) Review of the common fracture modes in laminates and their interaction 

(based on Lockheed and NASA studies) 

(2) Review of free-edge stress analysis and prediction of delamination sites 

(based on studies by Lockheed Palo Alto Research Laboratory and 
Drexel University) 

(3) Consideration of the problem of predicting delamination fracture in 

regions containing a high stress gradient (based on studies by Lockheed 
for NASA) 

(4) Analysis of delamination fracture by calculation of strain energy release 

rate (based on studies by Drexel University and Lockheed) 

(5) Consideration of the Interaction of interlaminar and intralaminar fracture 

(based on studies by Drexel University and Lockheed for the Air Force 
Office of Scientific Research) 

(6) Analysis of the interrelationship between strength and stiffness (based on 

studies by Lockheed for NASA) 

(7) Consideration of the correlation between strain energy analysis and 

statistical analysis approaches in predicting the criticality of a 
laminate crack (based on studies by Lockheed, Cornell University, and 
Drexel University) 

The following people contributed to the observations, analysis methods, and 
methodology of delamination fracture analysis described in this paper: J. G. 

Bjeletlch and W. J. Warren (Lockheed Palo Alto), A. S. D. Wang and G. E. Law 
(Drexel University), T. K. O'Brien (NASA Langley), and J. T. Ryder (Lockheed Rye 
Canyon) . Finally, many researchers have contributed to the understanding of 
delamination and its influence on laminate fracture. These people are acknowledged 
in the partial bibliography at the end of this paper. 
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Five cases that involve delamination fracture are presented. These are: 


(1) Free-edge delamination 

(2) Compression delamination 

(3) Delamination interaction with intralaminar cracking 

(a) Ply peeling 

(b) Delamination after transverse cracking 

(c) Simultaneous delamination/ transverse cracking 
These cases are illustrated in the five figures which follow. 


INTERLAMINAR FRACTURE - CASE 1 : FREE EDGE PLY DELAMINATION 




SINGULAR INTERLAMINAR STRESSES AT CRACK-TIP 
OPENING MODE SHEARING MODE 

ACTUAL CRACK IS IN A COMBINED MODE 
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INTERLAMINAR FRACTURE - CASE 2: COMPRESSIVE LAMINATE FAILURE MECHANISMS 



A 



CRACK TIP 


SINGULAR INTERLAMINAR CRACK TIP 
NORMAL STRESS 

CRACK GROWTH IS PRECIPITATED BY 
STRUCTURAL BUCKLING ACTION 

UNSTABLE CRACK GROWTH DETERMINES 
THE COMPRESIVE STRENGTH OF LAMINATE 


INTERLAMINAR FRACTURE - CASE 3 (a): PLY-PEELING (ALSO A CASE 

OF PLY DELAMINATION) 
(SURFACE DELAMINATION) 

A«x 



V 
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INTERLAMINAR FRACTURE - CASE 3(b): 

*FREE EDGE DELAMINATION (ALONG 25/90 PLANE. MIXED MODE) 



INTERLAMINAR FRACTURE - CASE 3(c): 

"DELAMINATION EMANATING FROM THE ROOT OF A TRANSVERSE CRACK (MIXED MODE) 
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ULTIMATE TENSILE STRENGTH VERSUS STACKING SEQUENCE 


Studies were conducted In the mid-1970's to examine the Influence of laminate 
stacking sequence on the tensile properties. The results shown below and in the next 
figure compare the ultimate strengths of a number of quasi-lsotropic laminates 
constructed of T300/934 tape and fabric (ref. 1). A large variation in strength was 
found. This variation was not predicted from laminated plate analysis where stacking 
sequence does not alter the in-plane ply stresses when the specimen is loaded under 
tension. 


TENSILE STRENGTH OF QUASI-ISOTROPIC T300/ 934 LAMINATES AS A FUNCTION 
OF PLY STACKING SEQUENCE (LAMINATES DESIGNATED HMF WERE FABRICATED 
FROM WOVEN GRAPHITE FABRIC) 

{From ref. 1} 


MATERIAL 

THICKNESS 

(IN.) 
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DELAMINATION AND TENSILE STRENGTH 


An examination of several independent studies confirms the observation that 
stacking sequences that are prone to delamination exhibit generally lower ultimate 
strengths. 


TENSILE STRENGTH VARIATION WITH STACKING SEQUENCE 

EIGHT-PLY, GRAPHITE-EPOXY. QUASI-ISOTROPIC LAMINATES 
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DELAMINATION 


NO DELAMINATION 
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WHAT CAUSES DELAMINATION 

Strain gages placed on tensile specimen edges to measure the "through-the- 
thlckness" normal strain (z-strain) as a function of tensile strain (x-straln) 
showed that the delamlnat Ion-prone laminates developed a negative Poisson ratio 
(e.g., the z-strain was positive during tensile loading). Thus interlaminar 
interaction that was unaccounted for in laminated plate analysis was providing a 
driving force for delamination. 


THROUGH-THICKNESS NORMAL STRAIN VERSUS APPLIED TENSILE STRAIN 
£x AS A FUNCTION OF T300/934 PANEL TYPE DESCRIBED IN PREVIOUS FIGURE 

iFrom ref. 1 } 
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WHERE DOES DELAMINATION OCCUR 


Examination of the edges of tensile specimens before and after failure revealed 
that the location of delamination depended on the particular stacking sequence and 
in some cases involved multiple delamination at several locations through the ply 
thickness. Transverse cracks were seen to initiate both before and after delami- 
nation, depending on stacking sequence. The extent of delamination was found to be 
greater at the edges than in the interior of the specimen. 


x-z EDGE VIEW OF (0/45/-45/90)^ LAMINATE AFTER LOADING TO 9S% UTS 

(From ref. 1 ) 
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EDGE VIEW OF (45/-45/0/90) LAMINATE AFTER LOADING TO 95% UTS 

(From ref. 1 ) 



LONGITUDINAL x-z SECTION OF (45/-45/0/90)^ TAKEN NEAR MIDPLANE 
BY SECTIONING AND POLISHING 

(From ref. 1 ) 
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EDGE VIEW OF (90/45/-45/0)^ LAMINATE AFTER LOADING TO 95% UTS 

(From ref. 1 ) 
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ANALYSIS OF INTERLAMINAR STRESSES 

Consideration of stress equilibrium and conservation of angular momentum led to 
an understanding of how interlaminar stresses arise and the development of simple 
models to approximate the level and sign of these stresses. 

ANISOTROPY IN STIFFNESS OR EXPANSION COEFFICIENTS 
CAUSES PLY OR LAYER STRESSES 
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IN-PLANE PLY STRESSES RESULT IN FREE-EDGE STRESSES 




202 




NUMERICAL ANALYSIS OF INTERLAMINAR STRESSES 



Finite-element and finite-difference analyses have been applied to models of an 
y— z section of a laminate under uniaxial loading in the x direction. The through- 
the-thickness Poisson ratio predicted from the analyses agreed well with experi- 
mentally measured values. 
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STACKING SEQUENCE DEPENDENCE OF THROUGH THICKNESS 
CONTRACTION IN T300/934 QUASI -ISOTROPIC 


LAMINATES (ref. 1) 




Panel No. 

Stacking 

Sequence 

V 

xz 

Experimental Calculated 

Ultimate 

Tensile 
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(ksi) 

IS 
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63 
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COMPARISON OF NUMERICAL METHODS 

The detailed stress analysis was found to depend somewhat on the fineness of 
the model discretization and the order of the assumed variation of displacement or 
stress field in the case of finite-element analysis. Stress gradients In both the 
y and z directions were found in these analyses, as shown In the figures that follow 


THIS STUDY 

PIPES AND PAGANO (ref. 2) 


FREE EDGE STRESSES IN [o, 90. 90. o]^ T 300/934 







PREDICTION OF DELAMINATION LOCATION 

As shown In the figures, the location of peak tensile stress In the analysis 
often agreed with the experimentally observed delamlnatlon location. 


DISTRIBUTION OF Oz NEAR FREE EDGE UNDER TENSILE LOAD 

IN X DIRECTION 


EIQ 





THE PROBLEM OF STRESS GRADIENTS 


A problem arose, however, when the numerical analysis of Interlaminar stresses 
associated with moisture gradients was considered. The figures show the moisture 
concentration as a function of y,z position at two different times during the 
moisture absorption process. At a time that was small with respect to the equilib- 
rium time, the gradient of concentration is very steep near the free edge. In the 
limit of very small times, this gradient is infinitely steep. Since absorbed 
moisture is known to swell the composite locally, a gradient of hygroscopically 
induced stresses is introduced. 

AVERAGE WEIGHT GAIN (LOSS) DURING ABSORPTION (DESORPTION) 

RELATIVE TO EQUILIBRIUM MOISTURE CONTENT M^ 

{From ref. 3} 



MOISTURE CONCENTRATION C/Cg AFTER ABSORPTION TIME 
(M/Ms«4%) ^ 


{From ref. 3f 
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PREDICTION OF DELAMINATION UNDER A STRESS GRADIENT 


The results of stress analysis In a unidirectional laminate during desorption of 
moisture from an initially saturated condition show that a transverse normal stress 
greater than the transverse tensile strength of the material is present during the 
first part of the desorption period, but gradually reduces to zero upon drying out. 
Although strength theories based on these local stresses predict delamination during 
desorption, no experimental evidence of delamination was found in a series of 
experiments conducted under a NASA Ames contract (ref. 3). It was recognized during 
the analysis of these results that stresses calculated at a distance from the free 
edge or interlaminar region which is on the order of the fiber diameter cannot be 
physically real. The analysis is based on the assumption that fiber and matrix 
phases are smeared together at the lamina level. That assumption is not appropriate 
on the microstructural scale. At this level, the interaction of individual fibers 
and intervening matrix must be considered. Thus the analysis of interlaminar 
stresses is inherently limited to a cross-sectional area which contains on the order 
of 10 fibers. In graphite-epoxy laminates this results in a minimum finite-element 
size of approximately 0.001 x 0.001 in. 

a, THROUGH-THICKNESS DISTRIBUTION NEAR 
FREE EDGE OF UNIDIRECTIONAL 8-PLY 
LAMINATE DURING DESORPTION FROM 
EQUILIBRIUM MOISTURE CONTENT 

{From ref. 

(MPa/M%) 


0 30 60 
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PLY THICKNESS AFFECTS DELAMINATION ONSET STRESS 


Rodlnl and Elsenmann (ref. 4) found that the onset stress In a series of 
quasl-lsotroplc laminates varied with ply thickness as approximately 1 over the 
square root of laminate thickness. Delamlnatlon stress analysis based on continuum 
theory Is Independent of the scale of the structure and was unable to predict the 
experimental results. Rodlnl et al. proposed that Weibull strength theory which Is 
dependent on specimen volume could be used to explain the results. 

A second approach, based on fracture mechanics, was taken at Drexel and Lockheed 
to account for the observed size dependence of delamlnatlon and transverse cracking 
onset in a series of tensile-loaded (25/-25/90n)s laminates where n ranged from 
1/2 to 8 (ref. 5). 


PLY THICKNESS EFFECTS 


T 300/5208 LAMINATES 

AT ONSET OF DELAMINATION 


EXPERIMENT 

PREDICTION 

(±45/0/90); 

45.3 ksl 



45.3 

49.0 ksi 


45.5 


(± 452 / 02 / 902 ) 

38.3 


s 

38.6 

36.0 


37.8 


(± 453 / 03 / 903 ) 

28.7 


s 

30.7 

29.0 
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CASE STUDY USING T300/934 GRAPHITE-EPOXY LAMINATES 


Laminate configuration: 

(±25/90n)s n = 4, 1. 2, 3. 4. 6. 8 

Loading condition: 

Uniaxial laminate tension .5^^ or 

Uniform curing temperature AT = T^ - T^ ^ 

o ■ » ^ 

Types of cracking studied: 

1. Transverse cracks in 90° layer 

2. Edge delamination along 90/90 interface (mode I) and 25/90 

interface (mixed mode), and transverse-crack-induced delamination 

3. Edge-delamination-induced transverse cracks 
Objectives of study: 

1. Onset and growth behavior of various modes of damage 

2. Their interrelations 

3. Physical and/or geometrical parameters influencing these behaviors 

4. Whether minimechanics approach is viable tool in studying damages in 

composite laminates 

5. How to bridge gap between micromechanics and macromechanics 









OBSERVED TENSILE STRAIN Ex FOR TRANSVERSE CRACKING, DELAMINATION, 
AND ULTIMATE FAILURE VERSUS NUMBER OF 90° PLIES IN T300/934 



12 3 4 


n 



ANALYSIS OF DELAMINATION STRAIN ENERGY RELEASE BY CRACK CLOSURE 


The crack closure method pioneered for delamination analysis by Rybicki et al. 
(ref. 6) was applied to the analysis of strain energy release rate associated with 
both delamination and transverse cracking. The analysis accounted for the con- 
tribution of thermal stresses to the total energy release rate. The crack closure 
method also allows the determination of the G values associated with each crack 
opening mode. 


Y.v 



Y,v 



■ a + Aa ' ■■ 

The work required to close the crack extension Is approximated by 

AW«»[F (u,-u ) + F (v -V ) + F (w -w )]/2 
xfg yfg zfg 

where F^,F^,F^ are the components of the nodal forces required to close nodes 
f and g together. Thus, the energy release rates for the three crack extension 
modes are approximated by: 

G_ F (v,-v )/2Aa 

I yfg 

hi - P^(Uf-“g>^2Aa 

F (w,-w )/2Aa 
III zfg 
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ENERGY ANALYSIS 


G = - ^ 

da Aa 

Energy released in crack opening equals work required to close crack: 
-Au = %F • D 

Let Fg, Dg be the elastic solution for a loading of = 10"® and let F 
be the elastic solution for a loading of AT = -1°F. 

Then 

® ' ilr + D^(iT)) 

This may be rewritten as 



THE CRITICAL FLAW CONCEPT 

The next figure below Introduces the concept of an effective flaw which 
preexists in a material prior to loading or one which grows subcritically in the 
sense of satisfying the Griffith critical stress criterion based on independently 
measured stiffness and toughness. 

If a critical flaw is not assumed to exist, then the analysis of transverse 
cracking of progressively thicker 90° layers in the (25/-25/90n)s laminate predicts 
that the energy release rate goes to infinity as the ply thickness is increased 
without limit. The consequence of this is that the predicted stress level for 
transverse cracking approaches zero for very thick 90° layers. This is inconsistent 
with the finite strength of transversely loaded unidirectional composites. 


METHOD OF LINEAR FRACTURE MECHANICS 

• CONCEPT OF ENERGY RELEASE RATE 

• G-CURVES AS FUNCTIONS OF CRACK EXTENSION 



FOR MODE-I COLLINEAR CRACKING IN ISOTROPIC ELASTIC SOLIDS. 
Gj = (tt/E) ao^ 

CRITERION FOR CRACK GROWTH: 



HENCE 

»c= y(<5cE)/K) 
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DETERMINATION OF THE CRITICAL FLAW SIZE 


To determine the critical flaw size a^, we postulate that It can be calculated 
by knowledge of the Independently measured elastic modulus E the fracture toughness 
6];^ and the unnotched tensile strength In the longitudinal and transverse direc- 
tions. For T300/934 composites, the substitution of these quantities Into the 
Griffith equation 


a 


c 



gives Sj, (transverse) = 0.010 In. and a^, (longitudinal) = 0.020 In. Since these 
are both much larger than the fiber diameter, our assumption of homogeneous ply 
properties In deriving the effective flaw size Is consistent with the absolute size 
of the flaw. The transverse critical flaw size was applied to the analysis of 
transverse and delamlnatlon crack Initiation in the analysis. The relative sizes 
of ac and am are important In assessing the stability of delamlnatlon growth, 
which is found to depend on the 90° ply thickness in the (25/-25/90n)s experiments. 

BEHAVIOR OF G-CURVES IN FREE-EDGE DELAMINATION 





CONSIDER l±em. 


EDGE DELAMINATION MAY OCCUR 


1. ALONG 90/90 INTERFACE 
(MODE I, OPENING) 

2. ALONG -0/90 INTERFACE 
(MODES I, m MIXED 



220 


FRACTURE MODE AND SEQUENCE PREDICTION 


Energy release rate analysis is qualitatively and quantitatively successful 
in predicting experimental behavior. 

OBSERVED TENSIII STRAIN Ex FOR TRANSVERSE CRACKING, DELAMINATION, 

AND ULTIMATE FAILURE VERSUS NUMBER OF 90O PLIES IN T30O/934 
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COMPARISON BETWEEN THEORY (DOTTED LINES) AND EXPERIMENT 
(PREDICTION BASED ON - 1.4; = 0.008" AND AT - 225°F) 



n 
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IMPROVING THE RESISTANCE TO DELAMINATION 


From the fracture mechanics analysis of delamination onset, it is clear that the 
onset stress for delamination can be raised if the critical energy release rate (or 
delamination toughness) of the material can be improved. The table given below 
shows (1) that order-of-magnltude Improvements in matrix toughness do not translate 
into large Improvements in the composite prepregged from tape, and (2) that 
composites fabricated from cloth prepreg show greater improvement in toughness 
than do the tape materials. The references presented emphasize the Importance of 
microstructural factors in controlling (and limiting) improvements in delamination 
toughness. 


FRACTURE ENERGIES OF GRAPH ITE -EPOXY MATRIX COMPOSITES 


Material 

Ref. 

Resin or 
Interlaminar 
Fracture 
Energy 
(KJ/m^) 

Comments 

5208 

7 

0.076 

TGMA/DDS Matrix 

Code 69/HTS Graphite Tape 

8 

0.076 

M 

5206/Morg II Graphite Tape 

9 

0.17 

II 

F263/T200 Graphite Cloth 

7 

0.36 

II 

5206/Morg 1 1 Tape 




UD Cross Fiber 

9 

20.4 


(+ 45)5 Cross Fiber 

9 

7.8 


Ql Cross Fiber 

9 

9.6 


A 

10 

0.16 

MY750 / NMA / BDMA 

A/HMS Tape 

10 

0.24 


B 

10 

0.33 

MY750 / Piperidine 

B/HMS Tape 

10 

0.28 


C 

10 

1.40 

MY750 / 3.2% CTBN 

C/HMS Tape 

10 

0.37 


D 

10 

2.2 

MY750 / 6.2% CTBN 

D/HMS Tape 

10 

0.36 


E 

10 

3.20 

MY750 / 9% CTBN 

E/HMS Tape 

10 

0.49 


205 

7 

0.27 

Diglycldyl ether bisphenol A/ 

205/T300 Graphite Cloth 

7 

0.60 

Novelac/Dicyandiamlne Catalyst 

205+CTBN 

7 

5.1 

13x12 3000 tow BOO cloth 

205+CTBN/T300 Graphite Cloth 

7 

4.6 

CTBN + HYCAR additives 
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STIFFNESS/STRENGTH RELATIONSHIPS 


A current Lockheed/NASA study is emphasizing the development of a structural 
analysis approach to assess the Influence of transverse cracking and delamlnatlon 
the resulting laminate stiffness and strength properties. 


OBJECTIVE 

• DEVELOP " STRUCTURAL” ANALYSIS FOR LAMINATES WITH 
ARBITRARY LAYUPS 

• PREDICT RELATIONSHIP OF MATRIX CRACKING AND DELAMINATION 
TO LAMINATE STIFFNESS AND TENSILE STRENGTH 

PREMISE 

• A LAMINATE IS A STRUCTURE 

• MATRIX CRACKS AND DELAMINATIONS MAY OR MAY NOT BE 
DAMAGE 

• A "STRUCTURAL” ANALYSIS OF ARBITRARY LAYUPS IS NEEDED 


SHEAR STIFFNESS REDUCTION 


Finite-element analysis was used to determine the effective in-plane transverse 
and shear moduli as a function of transverse crack density. The results indicate 
that for the experimentally determined crack densities, the loss of in-plane stiff- 
ness due to transverse cracking is very small. 


SHEAR STIFFNESS REDUCTION 
AS A FUNCTION OF CRACK SPACING 



S 
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LAMINA STIFFNESS LOSS AS A FUNCTION OF CRACK SPACING 

MASTER CURVES 



STIFFNESS LOSS DUE TO MATRIX PLY CRACKING 

T300 - 5208 GRAPHITE EPOXY 


LAYUP 

90° PLIES 
S^ ^22 

f22)o 

+ 45 O PLIES 

1 ®12 

' («12)o 

-45^' 

S 

h 

PLIES 

«12 

(<^ 12)0 

^ ^xx, % 

[0/90]j 

3.33 

65 

— 

— 

— 

— 

1.2 

[»± 45 ] j 

— 

— 

14.6 

95 

14.2 

95 

0.9 

[0/90/±45]^ 

2.37 

52 

13.6 

95 

— 

— 

2.2 
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STIFFNESS CHANGE WITH DELAMINATION GROWTH 

Correlation of delamlnatlon size and stiffness shows a much more significant 
effect on stiffness associated with delamlnatlon, which can be rather simply 
calculated from a modified laminate analysis. 

MEASUREMENT OF STIFFNESS CHANGE 
WITH DELAMINATION GROWTH 





RULE OF MIXTURES ANALYSIS OF STIFFNESS LOSS 
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STIFFNESS LOSS AND DELAMINATION SIZE 


A comparison between model prediction of stiffness loss and experimental data 
is given below. 


STIFFNESS LOSS AS A FUNCTION OF DELAMINATION SIZE 

[+45/0/90] T300-5208 GRAPHITE EPOXY 
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TRANSVERSE-CRACK EXTENSION 


Transverse cracking of a ply typically terminates at the interface of a ply of 
differing orientation. Microscopic and X-ray examination of the crack tip by 
researchers at VPI&SU and Lockheed, among others, has failed to reveal significant 
damage to the fibers in the neighboring ply, at least when the specimen has-been 
subjected to only static loading (private communication, K. L. Reifsnider). Recent 
analysis of the stress concentration and energy release rate associated with the tip 
of the transverse crack has confirmed the lack of damage potential associated with a 
transverse ply crack. 

The figure below shows the finite-element model discretization of a transverse- 
ply crack. The ply arrangement for the quasi-isotroplc configuration is shown, but 
three other layups were examined with the same element model. The stress concen- 
tration associated with a transverse 90°-ply crack is much less in the quasl- 
isotropic case than in the transverse (90) g laminate. This is a result of the 
relatively small load carried by the cracked ply in the quasi-isotropic laminate. 

The minimum element size used in this analysis contains approximately 10 fibers. 
Analysis of stresses with elements on a fiber scale would lose physical meaning, 
given the homogeneous-ply-stiffness assumptions of the analysis. 


LOCAL STRESS CONCENTRATION DUE TO MATRIX CRACK 

[0/90/ + 45 ]^ 




LAYUP 

SCF 

0/90/±45‘ 

s 

1.1 

0 / 0 / ± 45 ' 

s 

1.6 

[Oj, 

8 


1.7 

[ 90 ] 

8 


1.5 
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WEIBULL STRENGTH ANALYSIS 


Phoenix (ref. 11) has derived a model which can be used to obtain the mean 
strength of a unidirectional composite as a function of its voliune. The volume is 
given by the product of the number of fibers in cross section (n) and the length 
divided by the critical length (m) . For a typical unidirectional tensile specimen, 
nm is on the order of 10^. The figure below shows the characteristic strength 
distribution function of a single critical length of fiber and the shift of that 
function associated with a volume of nm = 10^. Taking W(x) = 0.5 as the mean 
strength of a given volume, one can construct a table of strength normalized to a 
single fiber strength or that of a typical tensile specimen. 


CHARACTERISTIC DISTRIBUTION FUNCTION Wlx) FOR BUNDLE STRENGTH 
AND ASSOCIATED WEIBULL APPROXIMATIONS (ref. 11). (ALSO SHOWN IS 
STRENGTH DISTRIBUTION FOR A BUNDLE WITH # FIBER ELEMENTS) 



0.2 0.3 0.4 0.5 0.6 0.7 


DIMENSIONLESS LOAD. x/Xg 
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TENSILE STRENGTH VERSUS VOLUME 


Given that the stress concentration factor for a transverse crack in a quasi- 
isotropic laminate is approximately 1,1 and the voliime associated with the finite 
element nearest the transverse crack tip is approximately 5 x 10^ in nm units, we 
can see from the volume dependence of strength derived in the previous table that 
"failure" of that element is not expected on the basis of the Phoenix theory 
(ref. 11). 


UNDIRECTIONAL TENSILE STRENGTH VS SPECIMEN VOLUME** 


CASE 


Ln (nxn) X/X(CASE A) X/X(CASE B) 


A, A single fiber of the 

critical length 0,0 

B, A tensile specimen of 

,080 X ,5 X 10, in, 17,9 

C, A single ply ,005 in, 
thick of length ,04 in, 
near tip of transverse 

crack 9,65 

D, A single layer of fibers 
of length ,04 in, at the 
tip of the transverse 

crack 7,13 


1,0 


,25 


,35 


.41 


4.0 


1.0 


1.4 


1.6 


** A fiber spacing of .0004 in. and a critical length of .04 in. 
were assumed. 
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STRAIN ENERGY RELEASE RATE ANALYSIS 



The figure presented below shows the results of strain energy release rate 
analysis at the tip of a transverse crack for the geometry shown in the inset. 
Several types of laminates were analyzed. The energy release rate is plotted as a 
function of tensile strain. The longitudinal fracture toughness of graphite-epoxy 
is on the order of 80 to 90 ft-lb/ln^. A transverse crack in a (02/902)s laminate 
does not release enough energy to allow the transverse crack to propagate across the 
0° ply. Thus the energy release rate analysis approach provides the same prediction 
of transverse crack stability as the Welbull strength analysis. 


STRAIN ENERGY RELEASE RATE ANALYSIS OF MATRIX CRACKS 
IN EIGHT PLY GRAPHITE EPOXY LAMINATES 
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FAILURE BY STRUCTURAL LOAD TRANSFER 

The analysis of the final fracture event for tens lie- loaded laminates does not 
appear to depend on the details of ply microcracking. As experimental observations^^ 
of laminate damage growth under tension have accumulated, the Importance of local 
delamlnatlons in transferring a local load increase to the major load-carrying plies 
had been postulated by Reif snider and colleagues (ref. 12). The schematic edge 
section of laminate damage under static (top) and fatigue (bottom) loading shows 
regions of the fatigued specimen where essentially all of the applied load must be 
routed into the 0° plies around regions of extensively delaminated and tranverse- 
cracked 90° and 45° plies. Prediction of these loads will require more detailed 
structural modeling. 


PREFRACTURE PAHERNS NEAR FRACTURE SITE (ref. 13) 
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"STRUCTURAL" FAILURE ANALYSIS 


The degree to which the details of crack interaction must be modeled to predict 
tensile failure are being investigated in several ongoing research programs. 
Structural analysis of composite response can be carried out routinely with numerical 
methods lite finite-element analysis at several different structural levels; 
mlcrostructural, laminate structural, and macrostructural. However, this does not 
provide the total solution ta fracture analysis of composites. Our understanding of 
the structural implications of cracking at each level must be Improved to provide 
the proper failure criteria for the structural analysis. 

FRACTURE OF STRUCTURAL COMPOS ITES 




Structural 

Mechanics 

Geometric 

imperfections 

Post-buckling 

Damage 


Macro- 

Mechanics 

Local Damage 
at Stress 
Concentrations 


Minl- 

Mechanlcs 

Ply Fracture 
Delamination 


Micro - 
Mechanics 

Fiber Fracture 
Matrix " 
Interface " 


FractureO==ONDE 


What is the significance of damage or defect on the structural integrity? 
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Him III iiiiiiii II I I 


I 

CONCLUSION AND DISCUSSION / 

1. Methods of classical fracture mechanics can be viable crack growth criteria 
for matrix-dominated sublamlnate cracks, such as transverse cracking and 
ply delamlnatlon. 

2. The energy release rate Increases (generally) with the size of crack, which Is 
constrained by the physical dimension of the laminate's ply structure, such as 
the thickness of the piles. Thus, ply thickness has a profound Influence on the 
growth process of sublamlnate cracking, and even on the final failure of the 
laminate. 

3. The f Inlte-element/crack-closure technique Is a useful computational tool to 
simulate sublamlnate crack growth, provided the crack paths are relatively 
simple. 

4. The DIB (dllodobutane) enhanced X-radlography Is a useful NDE tool for monitoring 
the growth of sublamlnate cracks. 

5. For T300/934, the value of Is In the range 1.2 to 1.5 In-lb/ln^, but other 

experiments found It In the range 0.5 to 0.8 (refs. 3 and 6). (See sketch.) 


~ 1 . 2 - 1. 5 


This difference should be verified by additional experimentation. 

6. For the best fit of data, the critical material crack size a^ is taken as 
0.008 in. This should be verified by micromechanical investigation. Generally, 
a^ is a function of the size of inhomogeneity (fiber), the stiffness property 
surrounding the crack, and other processing variables. However, it is not known 
whether a^ might also depend on the ply structure of the laminate. 

7. In the study of delamlnatlon, all predictions were based on the total energy 
release rate G = Gj + Gxi + Gju for both mode-I and mixed -mode cracks. 
Experimental correlation indicates a need for a mixed-mode crack criterion, 
such as 



^l/^lc 


^II^^IIc 


^IIl'^^IIIc 


= 1 


There Is no firm evidence that Gjjj. and actually exist, and experi- 

mental means to measure them Independently has not been developed. 
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Some of the observed cracks were actually *planar cracks'. In the case of a 
crack front with a complicated contour, the degree to which a 2-D crack growth 
criterion is readily applicable is not known. 
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